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Localized surface plasmons are confined collective oscillations of electrons in metallic nanopar¬ 
ticles. When driven by light, the optical response is dictated by geometrical parameters and the 
dielectric environment and plasmons are therefore extremely important for sensing applications. 
Plasmons in graphene disks have the additional beneht to be highly tunable via electrical stimula¬ 
tion. Mechanical vibrations create structural deformations in ways where the excitation of localized 
surface plasmons can be strongly modulated. We show that the spectral shift in such a scenario is 
determined by a complex interplay between the symmetry and shape of the modal vibrations and 
the plasmonic mode pattern. Tuning conhned modes of light in graphene via acoustic excitations, 
paves new avenues in shaping the sensitivity of plasmonic detectors, and in the enhancement of the 
interaction with optical emitters, such as molecules, for future nanophotonic devices. 


I. INTRODUCTION 

To concentrate an optical field in nanometric vol¬ 
umes is challenging but nanophotonic advances teach 
us how to control and collect light in the form of surface 
plasmons at such small scale. Collective plasmon os¬ 
cillations are confined at the interface between a metal 
and a dielectric and stem from the interaction of light 
and the conduction electrons in the metal. Localized 
surface plasmons (LSPs) on the other hand are coher¬ 
ent oscillations of the electric field generated when light 
waves are trapped within a conductive nanoparticle. 
The frequency of the electron oscillation does not only 
depend on the composition of the metal and its di¬ 
electric environment but the optical response is like¬ 
wise determined by the particle size and shape. This 
paradigm gives rise to ’’plasmonic colouring effects” as 
has been used for centuries in stained glass, which is 
being revisited recently,^ but it has also laid the foun¬ 
dations to highly efficient devices for biochemical and 
medical testing given the pronounced sensitivity of the 
LSP resonances. 

With the birth of graphene plasmonics new routes 
have been enabled to gain an even more profound con¬ 
trol for strong light-matter interactions. Graphene 
is known for its superior electrical,mechanical,^’^ 
thermal,^ and optical properties^’along with a poten¬ 
tial for novel applications and technology. Graphene 
plasmons hold many advantages as opposed to plas¬ 
mons in noble metals since they are more strongly 
confined.Graphene can be gated to constitute highly 
tunable and low-loss guides for light in future compact 
plasmon devices, which are useful for signal processing 
and infrared sensing.The linear dispersion of the 
Dirac fermions give rise to a broad landscape of tunable 
optical properties and collective plasmons oscillations. 
Hence, the implementation of graphene for tunable 
nanophotonic devices such as modulators, switches and 
photodetectors simply underpins the rapidly growing 


interest in this field. 

Apart from the geometrical layout in nanostructured 
graphene that determines the optical response, plas¬ 
mons in curved and nonplanar graphene geometries are 
also being considered.Likewise, electrically dop¬ 
ing the system with an applied electric field proves to 
be a versatile technique for tuning graphene nanodisk 
plasmons. Gontrolling LSPs and their mutual inter¬ 
actions via external stimulus, enabled control of the 
enhancement of the photon absorption with applica¬ 
tion to infrared opto-electric devices as recent exper¬ 
iments demonstrated.^^ Rather than tuning the opti¬ 
cal response via an electric field or chemical doping, 
current simulations show that mechanical vibrations 
in extended unpatterned graphene sheets create the 
necessary momentum to efficiently couple light into 
plasmons.These findings suggest that new plas¬ 
monic devices can be designed whose spectral sensitiv¬ 
ity, photoresponse, photon absorption and modulation 
speed could be tailored by sound. 

In this paper, we add a new ingredient to actively 
control localized graphene plasmons through mechan¬ 
ical deformations.^^'"^^ Instead of changing the prop¬ 
erties of LSPs through varying the charge carrier con¬ 
centration, deformation of graphene disks through vi¬ 
brations produces observable plasmon shifts. Further 
to this, we show that the acoustic modes can strongly 
modify the decay rates of nearby emitters,such as a 
quantum dot or excited molecules, when placed in the 
nearest proximity of a graphene disk as illustrated in 
Fig. 1 . In this context it is predicted that the modal 
shape of the vibrating nanodisk determines the un¬ 
precedented degree of the light-matter interaction en¬ 
hancement when stimulated by sound. The fabrica¬ 
tion of the resonator could be realized by depositing a 
graphene disk on a perforated Si02 layer and clamp¬ 
ing it to the edges by attaching a thin insulation layer 
with a hole on top of it. The acoustical excitation of 
graphene nanodisks can be practically realized by a fo- 
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cused laser with frequency modulated intensity. Puls¬ 
ing the laser generates phototermal excitations that 
create mechanical stress and cause the graphene nan¬ 
odisks to vibrate at the modulation frequency. 


II. MECHANICAL VIBRATIONS OF 
GRAPHENE NANODISKS 


representing conservation for both charge and energy. 
More specifically, p(r) is the induced surface charge 
density, $(r) and are the total and external 

potentials, a{uj) is the surface conductivity and, Cm is 
the averaged permitivity of the surrounding medium. 
For arbitrarily shaped graphene structures, Eq. (3) can 
be solved numerically as we are going to clarify in more 
detail below. 


A monolayer graphene nanodisk that is clamped at 
the outer rim can commonly be described as a pre¬ 
tensioned vibrating circular membrane.The gov¬ 
erning equation for time harmonic undamped trans¬ 
verse vibrations ^(,t) (along the 2 ; direction) reads 

= ( 1 ) 

Ps^ 

where /(O) is the external force stimulus of frequency 
Cl and, S the sur face a rea of the nanodisk. The speed 
of sound V = y^T/ps is very sensitive to the mem¬ 
brane preparation and can be altered by pre-stressing 
the nanodisk with a tension T by assuming an unal¬ 
tered surface mass density pg. In our investigation we 
choose an average tension of T = 0.2 N m“^ and the 
surface mass density ps = 2.96 x 10“^ kg m“^ as re¬ 
ported in the literature.In order to find the mechan¬ 
ical eigenmodes we seek the homogeneous solution to 
an undriven problem. Hence, with fixed edges this so¬ 
lution reads 

t) = AmnJm{kmnr) COs(to6>) sin(^^m„t), (2) 

where Amn is the complex modal amplitude, Jm the 
Bessel function of the first kind of order m, and kmn = 
3ran I^ whcrc 2ran Is the nth root of J^. From this ex¬ 
pression we get a linear relationship between frequency 
and disk radius, Ctmn = jmn'^lR- 


HI. OPTICAL EXCITATIONS OF GRAPHENE 
NANODISKS 


The optical response of graphene structures relies 
on treating the material as infinitely thin with a fi¬ 
nite conductivity. For the sufficiently large diameters 
we can safely ignore atomic details^^’^^ and nonlocal 
effects.^^”^^ In the electrostatic limit we seek the so¬ 
lutions of Maxwell’s equations (with time dependence 
through a self-consistent relation:^^’^^ 


f dr<^{r)p{r) = ^ f drV$(r) • V$ec.t{r) 

Js Js 

/ drV^(r) • V / dr‘ 

Js Js 
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IV. MODULATION OF LSP RESONANCES 
BY SOUND 


The electro-mechanical dynamics is a complex pro¬ 
cess, but thanks to the huge temporal mismatch be¬ 
tween the interacting resonances, we are able to sim¬ 
plify the problem. The wave oscillation of the LSPs (at 
THz) occur at much faster pace compared to the me¬ 
chanical modes (at GHz), which implies that the plas- 
mons will follow the moving membrane adiabatically. 
Thus, at a given time plasmons only notice a static 
deformation rather than a rapidly vibrating body. 

In practical terms, this means that we can initi¬ 
ate the numerical simulations by calculating the dis¬ 
placements of the acoustic modes, from which we get 
a static deformation of the graphene nanodisks that 
is needed to subsequently compute the electromag¬ 
netic (EM) response. In calculating the plasmonic 
properties, we solve Eq (3) through finite-elements, 
which concretely means that the potentials and in¬ 
duced charge densities are discretized over an acousti¬ 
cally deformed domain.In all simulations, we choose 
a Fermi level Ep = 0.4 eV (resembling typical doping 
levels achieved in experiments^^’^^’^^), and Sm = 1.355o 
(^0 vacuum permittivity) accounting for the screening 
effect from the surroundings but without any phonon 
contributions.^^ 

We begin with the assumption that the disk is ex¬ 
cited by a low intensity pulse such that the maximal 
disk displacement does not exceed the disk radius it¬ 
self. We work in a spectral range where three differ¬ 
ent modes can be excited; 4^01, ^ 02 , and 4^ii. Fixing 
the disk radius at R = 50 nm, these three modes vi¬ 
brate around f]^n/27r = 2.0 GHz, 4.6 GHz, and 3.2 
GHz respectively, corresponding to the following roots 
jf'oi — 2.4, jo 2 — 6.5, and jn cs 3.8. With these me¬ 
chanical resonance frequencies, we are on the safe side 
of our assumptions since the plasmon frequency of the 
nanodisks without an acoustical load is as high as 24 
THz. In Fig. 2(a) we excite the fundamental axisym- 
metric (0,1) mode with an amplitude of 50 nm. Since 
the mechanical cycle determines the EM interaction, 
we compute the plasmon frequency over the mechani¬ 
cal phase upto a quarter period = 7r/2. The elasto- 
static limit is met when ^ 0, meaning that the disks 
is not in motion and conventional LSPs in graphene 
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are excited. As seen in Fig. 2(a), we inspect solely 
the lowest breathing and dipole modes since higher or¬ 
dered resonances are very similar. When the mechan¬ 
ical phase accumulates, both resonances redshift but 
the breathing mode exhibits a more pronounced sensi¬ 
tivity as compared to the dipole resonance as depicted 
in Fig. 2a. The higher sensitivity results in a more no¬ 
table shift. Since the plasmonic response is modified 
by sound, the plasmo-mechanic coupling can be ap¬ 
preciated from the density of states (DOS) p{uj) of the 
plasmonic nanodisk while considering the acoustic load 
introduced by the mechanical vibrations. Note from 
Fig. 2a that p{uj) peaks around the band edges (tem¬ 
poral van Hove singularities), but also remains high 
within the band during the mechanical cycle. 

Next, we investigate the T 02 modal vibration, which 
again is axisymmetric but contain m = 2 nodal circles 
as illustrated by the inset of Fig. 2(b). As compared 
to the fundamental (0,1) resonance, the plasmonic 
breathing mode in the present example shows to ex¬ 
hibit a higher sensitive and is redshifted by almost 0.1 
eV upto the point where it degenerates with the dipole 
mode. We will try to shed some light on the spectral 
shift owing to the mechanical perturbation. With a rel¬ 
atively small displacement field T(r, t), we perform a 
Taylor expansion of 1 /1 r — r' + T (r, t) — T (r', t) | and can 
qualitatively give the energy shift S{Jiw) of the LSPs in¬ 
duced by the electro-mechanical coupling as follows 

S(Jiw) (X [ (irV^(r) • [ dr . VT(r, t). (4) 

Js Js |r-r'|2 

In Eq. (4), r and r' run along the x — y plane, and 
(3 = |T(r, t) — T(r', t)|/|r — r'|. The axiymmetric dis¬ 
placements of mode (0,1) and (0, 2) are seen to have 
only little influence on the plasmonic dipole mode. In 
these cases, the potential distributions and the induced 
charges, as visualized in the insets of Fig. 2, are mostly 
localized at the nanodisk edge when excited at the 
dipolar resonance. The EM energy is thus, predomi¬ 
nately concentrated along the mechanically stationary 
rim (|r| ^ R and \r'\ R) where the numerator of P 
goes to zero, causing the integral in Eq. (4) to be very 
small. In this context, as expected, the energy varia¬ 
tion 5{hLj) is more pronounced for the breathing mode. 
The charge distributes now throughout the bulk of the 
disks, which causes more efficient electro-mechanical 
coupling with more pronounced plasmonic sensitivity 
as measured by a large redshift. The order of the ax¬ 
isymmetric modes also play a significant role. The state 
To 2 has two nodal circles and Tqi only one. The dis¬ 
tance between positive and negative charges along z 
is thus enlarged for the (0, 2) mode due to the disk 
displacement and gives rise to a reduced Coulumb in¬ 
teraction and an increased redshift. 

We finally study the mechanical mode of the Tn 
state, which is the lowest non-axisymmetric vibration 


with the radial nodal line along the y axis as shown in 
the inset of Eig. 2(c). Interestingly, as we will see in the 
following, there is a certain interplay between this me¬ 
chanical polarization and the orientation of the dipolar 
mode. While the breathing modes behaves similar to 
the mechanical (0,1) mode from Eig. 2(a), two orthog¬ 
onal dipolar modes do now exist that are degenerate in 
the static limit Ot ^ 0. As the mechanical cycle sets 
in as depicted in Eig. 2(c), the dipolar mode along x 
exhibits a redshift while the ^—dipole mode blueshifts. 
The reduction of the Coulomb interaction in this case is 
also responsible for the redshift of the x—dipolar mode. 
The LSPs oscillate along the radial nodal line when the 
^—dipolar mode is excited. In this situation, as op¬ 
posed to the x-dipolar mode, the mechanical vibration 
does not increase the separation of opposite charges, 
but instead, it increases the Coulomb interaction so the 
energy increase causes the LSPs to blueshift. By chang¬ 
ing the frequency of the intensity modulated laser, we 
are able to excite many other mechanical modes, but 
we do not anticipate pronounced differences as com¬ 
pared to the previous discussion, which demonstrated 
an unique way to tune LSP resonances in graphene by 
acoustic stimulation. 


V. TUNABLE SPONTANEOUS EMISSION 

Modifying the photonic environment of emitters to 
enhance the spontaneous emission, has been realized 
by means of both dielectric structures^^’^^ and plas¬ 
monic cavities,in order to gain full control of the 
light-matter interaction. In what follows, we propose a 
novel strategy by imposing acoustic wave excitation to 
dictate the efficiency of the spontaneous emission en¬ 
hancement without altering the geometrical surround¬ 
ings. As we show in Eig. I, a dipole emitter is placed 
in the center at z = 60 nm in front of the nanodisk. 
Again, we concentrate on a set of specifically mechani¬ 
cally excited states and study their temporal evolution 
through the accumulated phase, as illustrated in Eig. 3. 
We begin by recalling the enhancement of the sponta¬ 
neous emission decay rate,^^ 

r _i , 67reo Im{//* ■ £;5(ro)} 

ro“^+ w ’ 

where k is the wave vector, p the dipole moment, and 
Es{ro) the scattering field at the dipole’s origin vq. The 
numerical results shown in Eig. 3 illustrate the complex 
interplay between the enhancement of the decay rate 
for specific dipole orientations and vibrational states. 
The peak locations of the Purcell factors are consis¬ 
tent with the eigenmodes calculated in Eig. 2, mean¬ 
ing that the resonances can be traced back to dipolar 
and breathing mode excitations of the nanodisks, when 
they are excited by an x—polarized and z—polarized 
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emitter, respectively. Overall, we can see that the en¬ 
hancement is remarkable with peak values approaching 
10 ^. 

In the following we will discuss on the selective en¬ 
hancement process with respect to orientation of the 
dipole and the mechanical cycle. The (0,1) mode 
makes the graphene disk move close to the emitter 
when the mechanical cycle flt/ir advances. As seen 
in Figs. 3(a) and 3(b), this is true for emitters, both 
polarized along the x and the 2 ; axis. At a higher acous¬ 
tic excitation frequency, the enhancement of the decay 
rate for a z—polarized emitter via the T 02 modal dis¬ 
placement, as seen in Fig. 3(d), can be traced back to 
the same mechanism as before. If however, a dipole res¬ 
onance is excited within the nanodisk where the electric 
charge confinement at the edge [see inset in Fig. 2(b)] 
moves away from the emitter due to the ( 0 , 2 ) modal 
shape, the Purcell factor decreases with time as illus¬ 
trated in Fig. 3(c). The mechanical vibration in the 
Til state breaks the axial symmetry and determines 
the optical response with dependence of the in-plane 
dipole orientation. For x—polarized emitters, the ex¬ 
cited dipole mode is effectively tilted by the mechanical 
displacement Tn. This effect, causes the edge of the 
nanodisk actually to move closer to the emitter with 
time and induces an enhancement of the decay rate as 
shown in Fig. 3(e). When polarized parallel to radial 
nodal line, the induced dipole remains stationary in po¬ 
sition, so the only effect is a slight blueshift as seen in 
Fig. 3(f). Finally, a complex interplay is predicted for 
the 2 ;—polarized case, giving rise to both Purcell factor 
enhancement and the excitation of x-polarized dipolar 
(the dotted elliptical regime) and breathing (the solid 
elliptical regime) resonances, see Fig. 3(g). 


VI. QUANTUM CALCULATIONS 

So far, our calculations of plasmons have been 
based on classical electrodynamics, Eq. (3), and mem¬ 
brane displacement entered only the plasmon dynam¬ 
ics through the modified Coulomb interactions associ¬ 
ated with the non-planar confinement of the electron 
gas. Here, we extend our study to include quantum 
considerations taking atomic-scale details and quan¬ 
tum plasmonics into account.We will consider the 
small-displacement limit, where the curved-space cor¬ 
rections to the Coulomb interaction can be neglected, 
while strain effects on the electron states are still ac¬ 
counted for. We start from a tight-binding (TB) rep¬ 
resentation of the flat membrane to first calculate elec¬ 
tronic states for the equilibrium case, i.e. in the ab¬ 
sence of electromagnetic fields. Subsequently, we in¬ 
voke the random-phase approximation (RPA) to cal¬ 
culate the energy-loss function associate with plasmon 
dynamics.For membrane dynamics we fill focus our 


attention on sufficiently small deformations where the 
in-plane lattice strains are captured by an empirical 
tight-binding parameter parametrizations,that work 
fairly accurately for a relative bond-stretching less than 
20%. In this model, the local tightbinding parameter 7 
reflects the local change in bond lengths a and follows 
the empirical expression^^ 

/^3.37a \ 

7 = 70 exp- 1 ( 6 ) 

V «o J 

with ao = 1.42Abeing the lattice constant of the un¬ 
strained lattice, while 70 = 2.8 eV.^^ 

For the purpose of illustration, we consider a disk of 
hexagonal shape with perfect six-fold rotational sym¬ 
metry and an armchair edge termination. As a con¬ 
sequence, the lattice supports no additional electronic 
edgestates at the Fermi level^^ and in accordance with 
group theory, both electronic states and plasmon states 
will be either nondegenerate or come has degenerate 
pairs. 

The hexagonal graphene disk that contains 366 car¬ 
bon atoms, being equivalent to a sidelength of roughly 
2nm. For the deflection amplitude Anm we limit 
ourselves to 0.4 nm, which corresponds to no carbon- 
carbond bonds being stretched more than 20 %. 

For our RPA calculations at room temperature 
{ksT = 25.7 meV) we follow previous work^^’^^ and 
use a plasmon damping of h/r = 1.6 meV and doping 
of the graphene disk corresponding to a Fermi energy 
of 0.4eV. 

In Fig. 4 we summarize the results which overall re¬ 
sembles the dynamics seen also from our classical elec¬ 
trodynamics considerations in Fig. 2. In particular, 
the plasmon breathing mode around 0.7 eV is signif¬ 
icantly modulated by the mechanical displacement of 
the membrane hosting the plasmons. In the example, 
none of the carbon-carbon bonds are stretched by more 
than 20 % (which is an upper critical limit for elastic 
deformation^^), yet the plasmon energy modulation is 
as high as 100 meV corresponding to a relative modu¬ 
lation of up to 40% for the dipole resonance while it is 
around 14% for the breathing mode. This suggest that 
plasmon-energy modulation can be appreciable. 

VII. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the modulation of 
LSPs in graphene nanodisks by mechanical vibrations, 
which is a new route to tune graphene plasmons. 

We note that classical electrodynamics have been 
used to illustrate the basic idea and for the ease of 
illustration, the membrane displacement is slightly ex¬ 
aggerated. However, even though our classical consid¬ 
erations will have their shortcomings, our full quan¬ 
tum treatment of plasmons in a realistically strained 
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tight-binding lattice confirm the existence of mechani¬ 
cal modulation of graphene plasmons. 

Rather than changing the shape and size of the disks 
or injecting charge carriers via doping or gating, ap¬ 
propriately chosen acoustic stimulations with respect 
to the modal profile, determines the very nature of the 
plasmonic response. This novel acousto-plasmonic sys¬ 
tem opens many striking possibilities to mold the fiow 
of light at the nanoscale and to enhance light-matter 
interactions. We have speculated how this can be used 
to modulated the light-matter interactions, such as the 
decay dynamics of a quantum emitter in the near vicin¬ 
ity of the membrane. In fact, such an exciting experi¬ 


ments has in the meantime been done.^"^ 
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dipole emitter 



X 2R 


PIG. 1 (color online). Schematic diagram of the system 
under investigation: A graphene nanodisk with radius R 
vibrates at its lowest asymmetric mode Tn. An illustrative 
example is depicted showing a quantum dot simulated as 
a x-polarized dipole emitter situated nearby the nanodisk. 
We compute the induced near-held potential as excited by 
the emitter. 
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Mechanical phase Qt / TT p(<^) Mechanical phase Qt / tt P((^) Mechanical phase Qt / tt P((^) 


FIG. 2 (color online). Plasmon frequency as a function of the mechanical phase O.t/'K for (a) ^oi, (b) ^ 02 , and (c) ^ 11 . 
The spectra show only the lowest energy plasmon resonances, which are the breathing (squares) and dipole modes (circles 
and triangles). The insets show the displacements in their respective mechanical states together with the potential 
distributions for the specific plasmon modes. 
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FIG. 3 (color online). Purcell factor for dipole emitters 
with X—, y— , and z— polarizations placed 60 nm away from 
the center (z = 0) of the vibrating nanodisk. The emission 
decay rate F is normalized to the vacuum rate Fo, and plot¬ 
ted against the photon emission energy at different values 
of the mechanical phase QZ/tt. We consider the following 
displacement fields: (a,b) Toi, (c,d) T 02 , and (e-g) Tn. 
The various dipole polarizations are indicated in the illus¬ 
tration. 
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FIG. 4 (color online). TB-RPA calculations of quantum 
plasmons in a 366 atoms graphene membrane with a me¬ 
chanical 4^02 displacement mode. Details of the parame¬ 
ters are discussed in the text. Note how the low-energy 
dipole modes are degenerate, while the high-energy breath¬ 
ing mode is nondegenerate. 




